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Technical  Report  Summary 

An  examination  of  the  spectral  characteristics  of  the  Lg 
phase  for  intra-plate  earthquakes  in  North  America  shows  a uni- 
form relationship  between  the  corner  period  and  the  seismic 
moment.  This  implies  a uniformity  of  earthquake  processes  in 
the  central  United  States  over  a wide  range  of  event  sizes.  Us- 
ing reduced  spectra  a method  has  been  outlined  for  relating  mag- 
nitude observations  at  periods  other  than  1 second  to  mb. 

Multi-mode  theory  is  applied  to  computing  ground  motion  time 
histories.  Agreement  between  predicted  and  observed  ground  motion 
shows  the  method  to  be  excellent  over  a distance  range  of  100  to 
2000  kilometers.  A comparison  of  predicted  ground  motions  com- 
puted for  a stable  continental  interior  model  and  for  an  active 
model  shows  that  effects  of  the  difference  in  the  attenuation  mod- 
els are  seen  at  distances  greater  than  100  kilometers  for  periods 
greater  than  1.5  seconds. 

Values  of  the  attenuation  coefficient  for  Rayleigh  waves  of 
period  5 to  40  seconds  have  been  determined  for  various  paths 
across  North  America  from  the  nuclear  events  Rulison  and  Rio  Blanco. 
Application  of  a Backus-Gilbert  inversion  method  showed  the  qual- 
ity factor  values  in  the  upper  crust  of  the  western  Cordillera  are 
about  half  those  of  the  upper  crust  in  eastern  North  America. 

Low  velocity  surface  sediments  strongly  affect  the  attenuation  at 
periods  less  than  5 seconds.  Similar  studies  applied  to  surface 
wave  traversing  Eurasian  paths  are  being  used  as  a tool  to  deter- 
mine attenuation  differences  by  tectonic  provinces  in  this  region. 
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Semi-Annual  Technical  Report 
Introduction 

The  investigations  under  study  in  this  project  pertain  to 
the  properties  of  sources  within  continental  interior  regions 
and  of  shallower  focus  thrust  faults  near  the  borders  of  plates 
as  these  affect  the  radiation  pattern  of  seismic  waves,  and  the 
attenuation  structure  of  continental  shields,  plateaus,  and 
younger  tectonic  provinces  along  typical  paths  from  such  sources. 
Both  areas  of  study  are  related  to  the  amplitude  of  seismic  waves 
and  magnitude  determination. 

Previous  Technical  Reports  have  presented  a study  of  the 
Ms  vs  mb  relation  for  central  U.  S.  earthquakes,  a technique  for 
determining  earthquake  source  parameters  by  using  surface  wave 
amplitude  spectra,  the  seismic  scaling  law  and  source  mechanism 
characteristics  for  small  magnitude  events  in  the  central  U.S. 
and  the  Q structure  above  the  subduction  zone  in  South  America. 

The  Report  for  the  present  period  will  concern  spectral 
characteristics  of  short  period  surface  waves  in  continental  in- 
teriors, the  Q structure  of  the  North  American  and  Eurasian  con- 
tinents, and  the  radiation  from  shallow  thrust  faults  in  eastern 
Siberia-Kurlle  Islands.  The  Report  will  include  both  work  that 
has  been  completed  and  that  in  progress. 
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A.  Work  Completed  During  the  Report  Period. 

Work  completed  during  this  six  month  report  period  has  been 
presented  in  four  papers  which  have  been  submitted  for  publica- 
tion . 

R.  L.  Street,  R.  B.  Herrmann,  and  0.  W.  Nuttli,  Spectral 
Characteristics  of  the  Lg  Wave  Generated  by  Central 
United  States  Earthquakes,  accepted  for  publication 
in  the  Geophysical  Journal . 

R.  B.  Herrmann  and  0.  W.  Nuctli,  Ground  Motion  ModelJing 
at  Regional  Distances  for  Earthquakes  in  a Continental 
Interior,  I.  Theory  and  Observations,  accepted  for  publi- 
cation in  Earth  Engineering  and  Structural  Dynamics. 

R.  B.  Herrmann  and  0.  W.  Nuttli,  Ground  Motion  Modelling 
at  Regional  Distances  for  Earthquakes  in  a Continental 
Interior,  II.  Effect  of  Focal  Depth,  Azimuth  and  Atten- 
uation, accepted  for  publication  in  Earth  Engineering 
and  Structural  Dynamics. 

B.  J.  Mitchell,  Regional  Rayleigh  Wave  Attenuation  in  North 
America,  submitted  to  J.  Geophys . Res. 

A digest,  or  presentation  of  aspects  of  these  papers  most  immedi- 
ately of  Interest  to  project  goals  follows. 

A . 1 . Spectral  Characteristics  of  the  Lg  Wave  Generated 
by  Central  U.  S.  Earthquakes, 
by  R.  L.  Street,  R.  B.  Herrmann,  and  0.  W.  Nuttli 

For  many  earthquakes  the  S phase  is  not  a distinct  arrival, 

i 


but  is  associated  with  other’  arrivals  of  body  and  surface  waves. 
The  theory  of  Brune  (1970,  1971)  and  of  Savage  (1972)  have  not 
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approached  the  problem  of  estimating  source  parameters  from  such 
a complicated  set  of  arrivals.  The  objective  of  this  study  Is 
to  accomplish  this  task. 

The  events  studied  occurred  in  the  central  U.  S.  between 
1961  and  mid  197^  and  fall  In  the  magnitude  range  0.5  < m^,  5.C. 

Sixty  three  seismograph  stations,  all  part  of  the  WWSSN,  LRSM, 
Saint  Louis  University,  or  independent  networks,  provided  Lg  ob- 
servations for  one  or  more  of  the  events  studied.  In  each  case 
the  instrument-corrected  spectra  of  the  ground  motion  were  used 
for  analysis  only  in  the  pass  band  for  which  the  instrument  re- 
sponse was  greater  than  one-tenth  the  maximum  response.  In  the 
digitizing  of  the  seismograms  the  resolution  attainable  depended 
on  the  time  rate  of  the  record.  LRSM  short  period  seismograms, 
when  viewed  through  a 10X  viewer,  were  found  to  give  good  spectral 
results  at  periods  as  short  as  0.1  seconds.  WWSSN  type  data  were 
useful  down  to  periods  of  0.3  seconds.  Records  were  generally 
digitized  from  the  time  of  the  maximum  trace  offset  associated 
with  the  Sg  or  Lg  arrival  until  the  signal  reached  the  background 
noise  level. 

Sample  displacement  spectral  densities  are  shown  in  Figures 
1 and  2.  For  both  events  the  spectrum  is  flat  at  the  long  period 
end  and  falls  off  at  very  short  periods  as  T^  (or  For 

some  events,  such  as  that  of  July  23,  1962  (Figure  2),  an  inter- 
mediate trend  was  observed 

The  shape  of  the  Lg  spectru-  in  Figures  1 and  2 is  similar 
to  that  predicted  by  3rune  (1970).  However,  some  different  in- 
terpretations of  the  spectra  are  required  since  body  wave  theory 
no  longer  is  applicable  over  the  propagation  paths  in  question 
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and  since  the  Lg  phase  is  not  a clear  S wave  but  most  likely  a 
superposition  of  higher  mode  surface  waves. 

A model  was  selected  which  properly  corrects  the  Lg  spectrum 
for  geometrical  spreading  and  at  the  same  time  relates  the  long- 
period  flat  portion  of  the  spectrum  to  the  seismic  moment  MQ.  The 
model  chosen  relates  the  Lg  spectrum  _fL(»o ) at  a distance  r from 
the  source  to  the  distance-corrected  far-field  estimate  of  the 
source  spectrum  S(C*J) 


f 4-rr(o^V0^Jra-'>  , V-SrD 

Uir 


Here  rQ  serves  two  functions.  First,  rQ  adjusts  the  scaling  so 
that  S(uj  =o)  = M0  . Secondly,  r0  establishes  a distance  which 
marks  the  change  in  the  character  of  the  geometrical  spreading 
of  the  Lg  wave  train  from  that  of  typically  body-wave  at  near  dis- 
tances to  typically  surface-wave  at  large  distances. 

The  values  adopted  for  the  parameters  in  equation  (1)  are 
P = 2.5  gni/cm^,  = 3.5  km/sec,  and  rQ  = 100  km.  The  value 
chosen  for  r0  is  that  which  was  found,  empirically,  to  give  stable 
values  of  M0  as  equation  (l)  is  evaluated  at  various  distances  r. 

Equation  (1)  was  applied  to  the  spectra  obtained  from  sev- 
enty-eight earthquakes  in  the  central  IJ.S.  For  each  event  average 
values  of  the  seismic  moment  and  corner  periods  were  determined 
from  the  distance-corrected  spectra . Figure  3 presents  the  gen- 
eral features  of  the  Lg  source  3pectra  estimates  for  representa- 
tive of  these  events.  It  is  noted  that  for  events  with  seismic 
moments  less  than  10^^  dyne  cm,  no  to  ^ trend  is  observed.  For 
events  v’ith  greater  3elsmic  moments,  (jo-*  trends  are  noted  in 


1 
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some  instances. 

Figure  4 is  a plot  of  the  To2  corner  period  versus  seismic 

O ~P 

moment,  where  To2  is  the  corner  period  between  the  eo  and  eo 
spectral  trends  (for  spectra  containing  both  a T0^  and  a T^2  cor_ 
ner,  the  mean  of  these  two  is  taken  as  Tq2),  The  average  trend 
of  the  corner  period  is  indicated  by  the  straight  line  segments. 
Special  features  of  this  plot  include  a)  the  consistent  relation- 
ship between  the  corner  period  and  the  seismic  moment,  Mq,  and 

b)  the  kink  in  the  cor-''*’  period  vs  Mq  curve  for  seismic  moments 
21  21 

between  10  x 10  and  5 x 10  dyne-cm. 

Several  conclusions  follow  from  the  features  that  are  noted 
in  the  plots  of  Figures  3 and  4. 

1.  The  consistent  relationship  between  TQ2  and  Mq  mark  a 
characteristic  of  earthquakes  of  the  central  U.S.  No 
such  relation  is  found  in  southern  California  where  wide 
ranges  in  stress  drop  introduce  a wide  variation  in  cor- 
ner period  for  a given  value  of  seismic  moment  (Thatcher 
and  Hanks,  1973). 

2.  The  kink  in  the  corner  period  is  a characteristic  of 
Aki's  (1972)  revised  model  A On  the  other  hand,  the  re- 
lation of  Mq  to  the  cube  of  the  corner  period,  in  regions 
excluding  the  kink,  is  a characteristic  of  Aki's  models 

A and  B,  but  not  his  revised  model  A. 

3.  Brune  (1970,  1971)  has  derived  a relation  between  the 
corner  period  To2  of  SH  waves  and  the  source  dimension 
whereby  under  conditions  of  constant  stress  drop  the 
seismic  moment  varies  as  the  cube  of  the  corner  period . 
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The  same  relation  applied  to  the  Lg  spectra  of  Figure  3 
imply  that  in  the  central  U.S.  earthquakes  with  M0  less 

pi 

than  10  dyne-cm  are  characterized  by  a stress  drop  of 
about  1 bar,  those  with  Mq  between  0 x 10^1  and  10^ 5 
dyne-cm  by  a stress  drop  of  about  6 bars  . The  kink  in 
Figure  U marks  a transition  between  the  two. 

4.  In  Figure  3,  in  terms  of  Brune's  (1970)  distinction  be- 
tween complete  and  partial  stress  drop,  in  the  central 

20 

U.S.  events  with  seismic  moment  less  than  3 * 10  dyne-cm 

_2 

exhibit  only  an  on  , or  complete  stress  drop,  whereas  some 
of  the  large  events  exhibit  also  an  or  partial  stress 

drop . 

Finally,  from  the  reduced  spectra  of  Figure  3 we  can  relate 
mp  (1  second)  to  magnitude  observations  at  periods  other  than  1 
second.  In  general,  the  equation  for  magnitude  for  a particular 
range  of  applicability  is  of  the  form  (Nuttli,  1972) 

nvj.  = B(T)  + C(T,D)  10g1Q(D)  + log10(A/T)  (2) 

where  D is  eplcentral  distance  and  A is  the  maximum  sustained 
ground  displacement  at  period  T of  a given  wave  type.  The  co- 
efficient C corrects  for  geometric  spreading  and  attenuation, 
the  factor  B scales  the  formula  to  source  reference  level. 
Assuming  that  the  magnitude  defined  for  a given  period  T is 
directly  related  to  log1Q  S(to),  a relationship  between  a mag- 
nitude scale  at  period  T and  some  other  period  T'  can  be  estab- 
lished by  comparing  the  relative  levels  of  S (cO)  at  the  two 
periods . 
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For  example,  from  Figure  3,  in  comparing  spectra  with 

po  PI 

M0  = 10  dyne-cm  to  MQ  = 10  dyne-cm,  we  would  expect  a 1.0 
unit  change  in  Ms  (20  second  period)  or  in  mb  (l  second)  but  a 
0.7  unit  change  in  magnitude  mg  ^ based  on  a 0.3  second  LG  com- 
ponent. Similarly  there  will  be  a 0.7  unit  change  in  mQ  ^ based 
on  a 0.1  second  component.  A scaling  law  can  thus  be  constructed 
giving  the  correspondence  between  M0,  Mg*  mb»  m0.3  and  mo.l 

Consequently,  in  any  region  such  as  the  central  U.S.  where 
source  characteristics  are  found  to  be  highly  predictable,  a re- 
liable estimate  of  mb  can  be  made  from  signals  with  periods  other 
than  1 second. 
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Figure  Captions 

Figure  1.  Spectrum  for  the  earthquake  of  02  February  1962, 
showing  the  T0^  corner  between  the  uJ3  and 
asymptotic  trends  of  the  spectrum. 

Figure  2.  Spectrum  for  the  earthquake  of  23  July  1962,  snow- 
ing the  T0i  corner  between  the  u>°  and  (ju-^  trends  and 

• _ 1 p 

the  T^2  corner  between  the  to  and  u>  trends. 


Figure  3 • General  shape  of  the  Lg  estimated  source  spectra  S(co). 
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A.  2.  Ground  Motion  Modelling  at  Regional  Distances 
for  Earthquakes  in  a Continental  Interior 

by  R.  B.  Herrmann  and  0.  W.  Nuttll 

I.  Theory  and  Observations 

This  study  attempts  to  model  the  ground  motion  at  regional 
distances  by  the  superposition  of  multiple  mode  surface  wave 
arrivals  . 

The  appi’ca-h  used  is  based  on  the  equivalence  of  a disloca- 
tion on  a fault  surface  within  an  elastic  isotropic  earth  model 
to  a double-couple  force  system  acting  at  the  point  of  application. 
This  equivalency  permits  one  to  calculate  the  theoretical  Fourier 
spectra  of  the  surface-wave  displacements,  using  normal  mode  sur- 
face wave  theory,  for  a dislocation  source  in  a multilayered  elas- 
tic isotropic  medium.  The  theoretical  development  and  the  result- 
ing equations  for  computing  the  ground  motion  are  given  by  Tsai 
and  Akl  (1970)  and  by  Levohln  and  yanson  (1971). 

In  addition  to  source  parameters,  theoretical  ground  motions 
depend  also  on  medium  effects,  such  as  larger  thicknesses,  com- 
pressional  and  shear  wave  velocities,  densities,  and  the  frequency- 
dependent  anelastic  attenuation  coefficients  of  both  fundamental 
and  higher  mode  surface  waves. 

Herrmann  has  shown  (1974 ) that  the  amplitude  spectrum  of 
the  surface  wave  ground  displacement  of  the  individual  modes  is 
relatively  insensitive  to  perturbations  in  the  velocities  or  den- 
sities of  the  layered  flat-earth  model.  However,  this  is  not 
true  of  the  phase  spectrum  of  the  ground  motion,  especially  at 
the  higher  frequencies.  As  a consequence  it  is  difficult  to 


lb 


obtain  a perfect  match  between  the  synthetic  seismograms,  for 
which  an  accurate  knowledge  of  the  phase  spectrum  is  required, 
and  the  actual  seismograms.  On  the  other  hand,  descriptive  prop- 
erties of  the  ground  motion  such  as  maximum  ground  velocity  and 
maximum  ground  acceleration,  which  depend  more  on  the  amplitude 
spectrum,  show  excellent  agreement  between  calculated  and  observed 
values,  as  will  be  demonstrated  below. 

Following  Tsai  and  Aki  (1970),  the  ground  motion  displace- 
ment of  the  surface  can  be  represented  in  the  form 

N(oj  ) 

f(r,t)  = 1 S(W)^_  A(r,w)  exp  (luJt  - ik  «r  - V'-r)  due  (l) 

7*  j-1  J I J 

In  this  equation,  S(uo)  is  the  Fourier  transform  of  the  disloca- 
tion time  history  and  Aj(r^uj)  exp(-  ikjr)  is  the  Fourier  transform 
of  the  medium  impulse  response  to  a point  dislocation  at  the  source. 
Aj(r,co)  is  a function  of  the  faulting  geometry  source  depth,  and 
the  elastic  properties  of  the  medium.  The  coefficient  takes 
into  account  anelastic  transmission  properties  of  the  medium. 

The  summation  runs  over  the  number  of  modes  present  at  a given 
frequency. 

The  numerical  computations  involved  in  determining  the 
Aj(r,eo)  at  high  frequencies  become  involved  and  time  consuming 
due  to  the  lack  of  numerical  accuracy  and  because  the  number  of 
modes  present  at  a given  frequency  is  proportional  to  the  fre- 
quency. In  the  earth  model  considered  below  we  were  limited  to 
a period  range  of  1.5  to  400  seconds  for  Love  waves,  which  had 
eight  inodes  present  at  a period  of  1.5  seconds,  and  to  a period 
range  of  2.0  to  400  seconds  for  Rayleigh  waves,  which  had  eight 
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modes  present  at  a period  of  2.0  seconds.  To  extend  the  compu- 
tations Into  the  period  range  of  0.5  to  1.5  or  2.0  seconds  in 
order  to  describe  the  known  large  surface-wave  motions  in  this 
period  range  would  have  required  the  use  of  about  thirty  sur- 
face-wave modes  at  a period  of  0.5  seconds. 


Computations 

In  applying  equation  (1)  a simplified  earth  model  for  the 
central  U.S.  based  on  Nuttli  et,  al_.  (1969),  modified  to  Include 
a low-velocity  layer  at  the  surface  and  a slight  low-velocity 
channel  in  the  crust,  wa3  used.  The  anelastlc  attenuation  model 
was  based  on  an  average  of  the  Rayleigh  and  Love  wave  attenuation 
values  determined  by  Mitchell  (1973a, b)  and  Herrmann  (1973,1974), 
and  on  the  higher-mode  short-period  Lg  values  by  Nuttli  (1973). 

Three  central  U.  S.  earthquakes,  whose  source  characteris- 


tics are  given  in  Table  1,  were  selected  to  test  the  appropri- 
ateness of  the  multiple-mode  surface  wave  theory  for  modelling 
ground  motions.  The  parameters  given  in  the  table  were  determined 
by  Herrmann  (1974)  and  by  Street  and  Herrmann  (1974). 

The  source  spectrum  model  is  based  on  the  models  of  Brune 
(1970)  and  Aki  (1972): 


S(tO ) 


Mo^i/iOo2)^ 


00 


(2) 


where  is  the  corner  frequency.  This  model  adequately  describes 
the  Fourier  amplitude  spectrum  of  the  source  time  function,  but 
may  be  in  error  with  respect  to  the  phase  information  at  hlgn  fre- 
quencies . 
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Comparison  of  Prodictions  and  Observations. 

In  this  study  we  are  only  modelling  surface  waves  with 
periods  greater  than  1.3  seconds.  Accordingly,  the  observational 
seismograms  which  we  shall  present  are  from  the  long  period  in- 
struments of  the  WWSSN  network.  These  instruments  peak  at  peri- 
ods of  10  to  30  seconds,  and  fall  off  at  6 db  per  octave  at 
periods  less  than  10  seconds  and  18  db  per  octave  for  periods 
greater  than  30  seconds.  Therefore  the  seismogram  embraces  the 
ground  displacements  in  the  period  range  of  10  to  30  seconds  at 
the  expense  of  the  motion  at  the  shorter  and  longe^  periods. 

Figure  1 compares  the  observed  and  predicted  long-period 
seismograms  for  the  vertical  component  Rayleigh  waves  at  FLO 
(Florissant,  Missouri)  and  the  Love  waves  at  OXF  (Oxford,  Missis- 
sippi) for  an  earthquake  which  occurred  near  Cairo,  Illinois  cn 
August  14,  1963.  The  distances  from  the  epicenter  to  FLO  and  OXF 
are  200  and  300  km,  respectively.  Note  the  excellent  agreement 
Detween  the  observed  and  predicted  seismograms  at  OXF  with  re- 
spect to  wave  form,  amplitudes,  and  signal  duration.  Likewise 
the  agreement  is  very  good  for  the  Rayleigh  wave  at  FLO.  Here 
the  important  motion  is  the  pulse-like  event  near  the  center  of 
the  record.  The  absence  of  the  shorter-period  tail  immediately 
following  the  large  pulse  on  the  predicted  seismogram  results 
from  too  high  a value  of  the  anelastic  attenuation  coefficient 
for  the  fundamental  mode  between  2 and  3 seconds. 

Figure  2 compares  the  observed  and  predicted  Love  waves 
at  DUG  (Dugway,  Utah),  at  a distance  of  2100  km  from  the  southern 
Illinois  earthquake  of  November  9,  1968.  The  similarity  in  the 
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large  amplitude  motion  is  quite  good,  although  the  predicted  Deak 
amplitudes  are  only  about  80%  of  the  observed.  The  latter  part 
of  the  observed  seismogram  is  not  accounted  for  by  pure  Love  waves, 
and  may  be  due  to  some  Rayleigh-wave  interference  or  multi  -path 
effects . 

Another  way  of  comparing  the  predicted  and  observed  ground 
motions  is  to  look  at  simple  diagnostic  factors,  such  as  maximum 
ground  velocity  or  acceleration  as  a function  of  distance.  For 
this  purpose  theoretical  velocitygrams  and  accelerograms  were  gen- 
erated for  various  distances  and  azimuths  about  the  source,  using 
the  appropriate  source  parameters  for  each  event  of  Table  l.  The 
results  are  given  in  the  full  paper. 

In  conclusion,  multiple-mode  surface  wave  theory  has  been 
used  effectively  to  describe  ground  motion  of  Love  waves  at  periods 
greater  than  1.5  seconds  and  of  Rayleigh  waves  at  periods  greater 
than  2.0  seconds.  We  feel  confident  of  being  able  to  predict 
ground  motions  at  distances  greater  than  100  km  from  the  sources 
on  the  basis  of  this  theory,  given  a suitable  model  of  the  earth- 
quake . 

II.  Effect  of  Focal  Depth,  Azimuth  and  Attenuation. 

In  further  development  of  the  ground  motion  studies  an  in- 
vestigation has  been  made  of  the  effects  of  focal  depth,  azimuth 
from  the  source,  and  of  attenuation  on  the  calculation  of  theo- 
retical ground  motion. 

Using  the  same  crustal  model  and  source  time  function  as  in 
Part  I,  two  anelaotic  attenuation  models  were  constructed:  a 
stable  model  appropriate  to  a geologically  stable  region  such  as 
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the  central  or  eastern  U.S.,  and  an  active  model  constructed  so 
as  to  be  characteristic  of  an  active  region  such  as  southern 
California.  At  shorter  periods  these  models  use  the  anelastic 
attenuation  coefficients  of  the  1-second  Lg  wave  estimated  by 
Nutt li  (1973)  for  the  central  U.S.  and  for  southern  California. 

For  periods  greater  than  10  seconds  both  models  are  assumed  equal 
and  are  based  on  Mitchell  (1973a, b)  and  Herrmann  (1973,1974),  as 
in  Part  I . 

For  the  source  model,  only  one  particular  fault  type  is  con- 
sidered, a pure  thrust  type  mot Jon  on  a fault  striking  north- 
south  and  dipping  45°  to  the  east  or  west.  From  Tsai  and  Aki 
(1970),  for  this  source  configuration  the  Rayleigh  spectral  ampli- 
tude will  vary  with  azimuth  ^ as 

A + B cos  2 ^ (1 ) 

where  A and  B are  functions  dependent  on  focal  depth,  mode,  and 
frequency.  The  radiation  pattern  at  a given  period  and  mode  may 
be  circular,  elliptical,  or  a two  or  four-leaved  rose,  or  com- 
bination of  these.  The  azimuthal  dependence  of  the  Love  wave  is 

C sin  2 (2) 

where  C is  a similar  function  dependent  on  depth  mode  or  fre- 
quency. Here  the  radiation  pattern  is  always  a four-leaved  rose. 

The  effect  of  focal  depth  on  the  amplitude  spectrum  of  Love 
waves  at  a distance  of  1000  km  along  an  azimuth  of  49°  is  shown 
in  Figure  3.  The  solid  curve  is  the  amplitude  spectrum  for  a 
source  depth  of  2.5  km,  the  dashed  curve  for  19*5  km.  But  the 
difference  in  focal  depth  causes  a profound  change  in  the  Love 
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wave  excitation  at  short  periods.  The  shallower  event  will  ex- 
cite the  shorter  periods  about  five  times  more  than  the  deeper 
event . 

Figure  4 shows  the  effect  of  focal  depth  for  Rayleigh  waves. 
The  spectral  hole  at  6 second  period  for  h - 2.5  km*  and  at  40 
seconds  for  h = 19-5  are  characteristic  of  the  choice  of  fault 
type  modeled  (Tsai  and  Aki,  1970).  Similarly  as  for  Love  waves, 
the  short  period  excitation  is  about  five  times  greater  for  the 
shallow  event  than  for  the  deeper. 

The  effect  of  the  stable  and  active  attenuation  models  on 
the  ground  displacements  of  the  Love  wave  at  distances  of  50  to 
500  km  along  an  azimuth  of  45°  for  a source  at  depth  9-5  km  is 
shown  in  Figure  5.  Differences  arising  from  the  attenuation  model 
are  not  measurable  at  distances  less  than  100  km.  This  is  for  the 
periods  2 seconds  and  greater  of  our  cdiuclatlons . At  large  dis- 
tances the  shorter  periods  are  attenuated  notably  more  by  the 
active  model  than  by  the  stable . 

Figure  6 compares  the  effects  of  the  two  attenuation  models 
on  the  vertical  component  Rayleigh  wave  ground  displacements. 

The  Rayleigh  displacements  are  not  as  impulsive  as  the  Love  wave 
values.  But  the  effects  of  the  different  attenuation  models  are 
similar  in  the  sense  that  the  active  model  filters  out  short 
period  components  at  distances  greater  than  100  km. 

Ground  velocities  and  accelerations  for  both  Love  and  Ray- 
leigh wave  motion  have  also  been  computed  and  compared  for  the 
two  attenuation  models. 

The  effect  focal  depth  is  further  examined  for  Love  and 
Rayleigh  wave  displacements  and  velocities  in  Figures  7 and  6, 


respectively.  The  motions  shown  were  generated  at  a distance  of 
250  km  from  the  source  along  a 1 azimuth  of  45°,  and  for  focal 
depths  of  2.5,  9.5,  and  19.5  km.  As  might  be  expected  from  the 
spectra  of  Figures  3 and  4,  the  short  period  excitation  Is 
strongly  dependent  on  focal  depth.  Note  also  that  a shallow 
event  has  a longer  signal  duration  than  a deeper  event,  due 
primarily  to  a greater  excitation  of  the  fundamental  mode  waves 
at  3hort  periods  which  in  turn  have  lower  group  velocities  than 
the  higher  mode  waves  at  the  same  period  (Panza  e^  a_L.,  1973)- 
For  the  particular  source  geometry  4 n question,  no  notable 
effect  of  azimuth,  other  than  variation  in  relative  amplitude. 

Is  found . 
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Figure  Captions 

Figure  1.  Comparison  of  observed  and  predicted  waveforms  re- 
corded on  the  OXF  and  FLO  long-period  instruments  for  the 
earthquakes  of  14  August  1965. 

Figure  2.  Comparison  of  the  observed  and  predicted  Love  waves 
on  the  DuG  long-period  seismograph  for  the  earthquake  of 
9 November  1966. 

Figure  3-  Love  wave  amplitude  spectra  for  focal  depths  of  2.5  km 
(solid  curve)  and  19*5  km  (dashed  curve). 

Figure  4.  Rayleigh  wave  amplitude  spectra  for  focal  depths  of 

2.5  km  (solid  curve)  and  19.5  km  (dashed  curve). 

Figure  5-  Love  wave  displacement  time  histories  as  a function 

of  distance  and  attenuation  model  for  a focal  depth  of  9.5  km 
and  receiver  azimuth  of  45°.  Displacements  are  in  cm. 

Figure  6.  Rayleigh  wave  displacement  time  histories  as  a func- 
tion of  distance  and  attenuation  model  for  a focal  depth  of 

9.5  km  and  a receiver  azimuth  of  45°. 

Figure  7-  Effect  of  focal  depth  variation  on  Love  wave  displace- 
ments and  velocities  for  the  stable  attenuation  model  at  a 
distance  of  250  km  and  an  azimuth  of  45°  from  the  source. 

Figure  8.  Effect  of  focal  depth  variation  on  Rayleigh  wave 
displacements  and  velocities  for  the  stable  attenuation 
model  at  a distance  of  250  km  and  at  an  azimuth  of  45°  from 
the  source  . 
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A.  3-  Regional  Rayleigh  Wave  Attenuation  in 
North  America 
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by  Brian  J.  Mitchell 
Introduction 

Surface  wave  velocities  and  their  regional  variations  have 
been  studied  extensively.  Pew  studies  have  been  made,  however,  of 
the  attenuation  of  these  waves,  and  there  is  almost  no  information 
on  the  regional  variation  of  surface  wave  attenuation. 

The  present  study  seeks  to  investigate  regional  variations 
in  surface  wave  attenuation  by  considering  eastern  and  western 
North  America  separately.  Various  studies,  including  measurements 
of  travel  times  (Hales  et  al,  1968),  geomagnetic  variations  (Por- 
ath,  1971),  and  heat  flow  (Roy  et  al,  1972)  indicate  a rapid  lat- 
eral change  in  properties  of  the  crust  and  upper  mantle  near  a 
line  separating  the  Cordillera  and  Great  Plains  of  North  America. 
Two  nuclear  events,  Rulison  and  Rio  Blanco,  were  detonated  in 
western  Colorado,  not  greatly  distant  from  that  boundary.  Ray- 
leigh waves  propagating  eastward  and  westward  from  those  events 
provide  a means  for  Investigating  differences  in  Rayleigh  wave 
attenuation  for  the  two  provinces  described  above. 

Effect  of  Velocity  Structure  on  Attenuation. 

Anderson  et  al  (1965)  present  a theory  for  the  attenuation 
of  dispersed  wavj  . They  determine  quality  factors  for  Love  and 
Rayleigh  waves  in  terms  of  the  quality  factors  for  compresslonal 
and  shear  waves  and  partial  derivatives  of  Love  and  Rayleigh  wave 
phase  velocities  with  respect  to  shear  and  compresslonal  wave 
velocities.  Surface  wave  attenuation  coefficient  values  are 
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related  to  the  surface  wave  Q values  through  the  equation 
y = if/UTQ,  where  y Is  the  surface  wave  attenuation  coefficient, 

U Is  the  surface  wave  group  velocity,  Q Is  the  surface  wave  qual- 
ity factor,  and  T is  the  period.  y , U,  and  Q may  refer  to  either 
Love  or  Rayleigh  waves. 

It  is  convenient  for  the  present  problem  to  reformulate  the 
equations  of  Anderson  et  al  (1965)  in  terms  of  surface  wave  atten- 
uation coefficients  rather  than  quality  factors.  The  equations 
then  become 


where  and  ^ are  compresslonal-  and  shear-wave  velocities, 
and  are  the  quality  factors  for  compresslonal  and  shear 

waves,  CL  and  Cp  are  Love  and  Rayleigh  wave  phase  velocities,  ana 
Si  is  a layer  index.  The  subscripts  to  , CC  , , and  refer  to 

the  frequency,  compresslonal  velocity,  shear  velocity,  or  density 
being  held  constant. 

Partial  derivatives  of  Rayleigh  wave  phase  velocity  with  re- 
spect to  compresslonal  and  shear  velocities  were  computed  for 
models  taken  to  be  representative  of  the  central  United  States 
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(McEvilly,  1965)  and  the  Colorado  Plateau  (Bucher  and  Smith, 

1971).  Figure  3 depicts  these  two  models  and  the  Rayleigh  wave 
phase  velocities  which  they  produce.  The  greatest  differences 
in  the  phase  velocities  occur  at  periods  less  than  5 seconds,  the 
Colorado  Plateau  model  being  much  slower  than  the  central  United 
States  model  for  those  periods.  At  greater  periods,  the  Colorado 
Plateau  velocities  parallel  and  lie  somewhat  below  those  of  the 
central  United  States . 

The  application  of  equations  (l)  and  (2)  requires  the  par- 
tial derivatives  discussed  above  and  appropriate  values  for  Q ^ 
and  Qp  . Two  simple  2-layer  Q models  were  selected  and  Rayleigh 
wave  attenuation  coefficients  were  computed  using  equation  (2). 

At  depths  greater  than  17  km  both  models  possessed  constant  Qp 
values  of  2000.  At  shallow  depths,  one  model  consisted  of  mater- 
ial with  a value  of  200  and  the  other  had  a Qp  value  of  100. 
These  are  simplified  forms  of  an  earlier  model  determined  by 
Mitchell  (1973^).  The  values  were  taken  to  be  twice  as  great 
as  the  values,  using  the  relationships  determined  by  Anderson 
et  al  (1963)  for  the  mantle,  and  assuming  they  also  apply  to  the 
crust.  Surface  waves  are  much  less  sensitive  to  than  to  . 

Knowledge  of  the  attenuation  coefficients  permits  one  to 
calculate  the  effect  on  amplitudes,  since  Ac*.  L“  )f(u  • 

The  distance,  x,  was  taken  to  be  1000  km,  and  the  amplitude  fac- 
tor, » was  computed  for  the  Colorado  Plateau  and 

central  United  States  velocity  models  (Figure  3),  combined  with 
the  two  Q models  described  above . Figure  4 presents  the  results 
of  the  computations.  Although  the  velocity  structure  produces 
differences  of  more  than  2 orders  of  magnitude  at  a period  of  1 
second,  it  has  very  little  effect  at  periods  of  3 seconds  and  more 
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Although  the  two  velocity  models  differ  throughout  their 


depth,  the  greatest  differences  occur  In  the  uppermost  few  km, 


where  low  velocity  sediments  are  present  In  the  Colorado  Plateau 


model  and  absent  in  the  central  United  Stater,  model.  The  short- 


period  Rayleigh  wave  energy  travelling  across  the  Colorado  Pla- 


teau model  concentrates  in  the  faster,  deeper  material  rather  than 


the  low-velocity  sediments 


The  Rayleigh  wave  radiation  patterns  at  a period  of  2 seconds 


In  Figures  1 and  2 indicate  very  large  amplitudes  to  the  north 


and  northeast  and  smaller  amplitudes  to  the  west,  south  and  south- 


east. The  large  amplitudes  were  recorded  at  stations  In  the  Canad- 


ian shield  and  northeastern  United  States  where  sediments  are 


either  thin  or  absent,  whereas  smaller  amplitudes  were  recorded 


by  stations  with  thicker  accumulations  of  sediments,  such  as  the 


western  United  States  and  Gulf  coastal  plain. 


Rayleigh  Wave  Attenuation  Coefficients 


For  a near  surface  composite  of  an  explosion  and  double-couple 


source  with  approximately  the  same  source-time  functions,  the  far- 


fleld  vertical  displacement  for  a Rayleigh  wave  can  be  written 


U;  (tu\ 


where  We  represents  the  vertical  displacement  produced  by  an  ex- 


plosion, F Is  the  relative  strength  of  the  double  couple,  ©G  Is 


the  orientation  of  the  double  couple,  ^ Is  the  Rayleigh  wave 
attenuation  coefficient,  and  x and  © are  the  distance  and  azimuth 


from  the  source,  respectively  (Toksoz  et  al,  1971). 


Values  for  We,  F,  ©Q  and  ^ can  determined  at  each  period 
by  fitting  equation  (3)  to  the  observed  data  by  a least-squares 
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iterative  process.  Equation  (3)  strictly  refers  to  a composite 

! 

of  an  explosion  and  a horizontal  double-couple  source.  The  hori- 
zontal double-couple  can  be  taken  to  represent  a vertical  strike- 
slip  fault  (Toksoz  and  Kehrer,  1972),  in  which  case  0O  will  oe 
the  strike  of  the  fault  (or  the  auxiliary  plane)  and  P will  be 
the  amount  of  tectonic  strain  release  relative  to  the  explosion. 

Although  equation  (3)  represents  a somewhat  special  case,  the 
amplitude  radiation  patterns  can  assume  a variety  of  shapes,  de- 
pending on  the  values  of  F and  0O  . The  patterns  are  elliptical 
or  two-lobed  when  F is  less  than  1.0,  and  four-lobed  when  F is 
greater  than  1.0.  Considering  the  scatter  in  amplitude  data, 
most  theoretical  and  observed  patterns  can  be  brought  into  satis- 
factory agreement  with  appropriate  values  of  F and  0O  . The 
attenuation  coefficient  value,  f , can  consequently  be  correct 
even  though  the  representation  of  the  source  as  a composite  ex- 
plosion and  vertical  strike-slip  fault  may  be  incorrect. 

The  mean  determinations  of  ^ , We,  ^ anc*  ^o  anc*  their  stan- 
dard deviations  for  data  from  Rulison  and  Rio  Blanco  appear  in 
Figures  5 and  6.  The  fit  of  the  theoretical  patterns  to  the  ob- 
servations, using  these  values,  are  shown  in  Figures  1 and  2. 

We  are  most  Interested  in  values  of  the  Rayleigh  wave  atten- 

-3  - 1 

uation  coefficient,  y . They  decrease  from  about  O.b  x 10  *m 
at  a period  of  5 seconds  to  very  small  or  slightly  negative  val- 
ues at  periods  of  12-14  seconds,  and  assume  a roughly  constant 
value  of  about  0.0b  x 10"3  at  greater  periods.  Results  of  the 
following  section  will  indicate  that  regional  variations  in 
structure  substantially  affect  Rayleigh  wave  attenuation  to 
periods  at  least  as  great  as  l6  seconds.  The  regional  variations 
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are  reflected  by  large  standard  deviations,  however,  only  at 
periods  less  than  fa  or  9 seconds  . Non-random  as  well  as  random 
effects  must  therefore  be  operative  on  the  observations  at  periods 
less  than  about  lb  seconds . Non-random  effects  probably  produce 
the  negative  attenuation  coefficient  values  at  periods  of  12-14 
seconds.  This  can  be  explained  in  terms  of  the  average  path 
lengths  from  the  sources  to  eastern  or  western  North  American 

B 

stations.  On  the  average,  paths  to  western  stations  are  substan- 
* tially  shorter  than  paths  to  eastern  stations.  Since  attenuation 

is  greater  to  the  west  than  to  the  east  from  the  sources,  tne  net 
effect  is  to  produce  low  apparent  attenuation  coefficient  values. 
Consequently,  all  of  the  observed  attenuation  coefficient  values 
for  periods  less  than  16  seconds  will  be  lower  than  the  actual 
average  values. 

Comparative  Attenuation  for  Eastern  and 
Western  North  America 

It  was  described  earlier  how  we  would  divide  North  America 
into  two  provinces,  one  east  of  the  Cordillera-Great  Plains  bound- 
ary, and  the  other  to  the  west.  The  amplitudes  recorded  for  sur- 
face waves  which  have  travelled  entirely  within  each  of  these 
provinces  are  now  considered  separately.  Disregarding  those  ob- 
servations which  do  not  clearly  represent  one  of  the  two  prov- 
inces, average  amplitudes  are  determined  for  the  eastern  (A^) 
and  western  province  (AgJ  for  each  period.  The  average  ratios 
of  the  values,  A-^Ag,  appear  in  Figure  7.  At  periods  of  lfa  sec- 
onds and  greater,  the  values  are  only  slightly  less  than  unity, 
at  shorter  periods,  the  values  decrease  with  decreasing  period 
to  a value  of  about  0.2  at  a period  of  7 seconds. 
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These  amplitude  ratios  lead  to  values  for  attenuation  coef- 
ficient differences  through  the  relation 


yc 


(4) 


where  and  Ye  are  average  attenuation  coefficient  values  for 
western  and  eastern  North  America,  respectively.  Differences  in 
anelasticity  models  can  then  be  investigated  by  equations 
analogous  to  (3)  and  (4).  For  differences  in  Rayleigh  wave  atten- 
uation coefficients  we  can  write 


(5) 


« 


The  validity  of  this  approximation  requires  that  the  first  fac- 
tors in  each  term  have  about  the  same  values  for  the  two  models 
'onsldered.  For  periods  of  5 seconds  and  greater,  equation  (5) 


was  lound  to  be  applicable.  If  we  assume  that  Q • * 2 Q 


discussed  above,  equation  (5)  becomes 


t ’ 


as 


This  equation  permits  a determination  of  the  difference  in  the 
vp,  distribution  for  two  different  models  using  the  data  of  Fig- 
ure 7.  As  seated  above,  we  will  consider  the  two  models  to  be 
the  two  large  provinces  of  North  America,  overlooking  smaller 
scale  variations  within  each  province. 
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The  basic  theory  for  Inverse  problems  permitting  the  consid- 
eration of  errors  and  resolution  was  developed  by  Backus  and  oil- 
bert  (1970).  The  present  Inversion  problem  is  particularly  simple 
since  it  is  a linear  one.  The  method  of  solution  is  discussed 


by  Mitchell  (1973b). 
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The  results  of  two  inversions  appear  in  Figure  8.  The  solid 
line  permits  a very  good  fit  to  the  data  of  Figure  7>  whereas  the 
dashed  line  is  a solution  for  which  the  theoretical  values  differ 
from  the  oDserved  values  by  about  30  per  cent  at  short  periods. 

The  main  features  of  both  curves  are  the  same . The  resolving  ker- 
nels on  the  right  side  of  Figure  ti  indicate  the  degree  of  resolu- 
tion for  A (Qp,1)  at  the  indicated  depth  in  the  model.  Narrow 
peaks,  well-localized  at  the  proper  depths,  indicated  by  the  ver- 
tical dashes,  imply  a high  degree  of  resolution,  whereas  broad, 
poorly-localized  peaks  imply  poor  resolution.  At  a depth  of  2 km, 
the  peak  is  narrow,  but  poorly  localized.  The  localization  im- 
proves below  2 km,  but  the  kernels  broaden  rapidly  at  depths 
greater  than  15  km. 

If  a Qp  model  for  eastern  North  America  is  available,  or 
assumed,  one  can  determine  an  average  model  for  western  North 
America  by  using  the  results  presented  in  Figure  8.  The  only 


crustal  model  available  at  the  present  time  is  that  of  Mitch- 
ell (1973b).  The  solid  line  in  Figure  9 is  a slight  simplifica- 
tion of  that  model  for  eastern  North  America.  The  dashed  line  in 
Figure  9 Is  a model  for  western  North  America  determined  from 
a modification  of  the  ENA  7 model  of  Mitchell  (1973b)  and  the 
values  indicated  in  Figure  8.  The  average  values  of  A or  the 
upper  crust  of  western  North  America  are  about  half  those  of 
model  ENA  7-  At  g’ eater  depths,  the  model  is  less  reliable 
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because  of  the  poor  resolution  in  the  determination  of  ^ 
at  those  depths. 

Effect  of  Temperature  on  Attenuation. 

One  of  the  results  of  the  present  study  is  that  values  for 
in  the  upper  crust  of  North  America,  in  gross  terms,  are  about 
half  those  of  eastern  North  America.  The  reason  for  this  differ- 
ence is  not  clearly  known,  although  temperature  differences  for 
the  two  regions  provide  one  possible  mechanism.  This  possibility 
can  be  examined  by  considering  the  variation  of  Qp  in  a homo- 
geneous material  a3  being  due  to  variations  of  temperature  (T) 
and  pressure  (P)  only.  This  situation  can  be  expressed  by 

qj'O-.r)  - <?;'  e>P[E'(r-r„VRtrl^p  [-PO'er]  (?) 

where  E*  and  V*  are  the  activation  energy  and  activation  volume, 
and  Qq1  , and  TQ  refer  to  reference  values  of  quality  factor  and 
temperature.  If  we  consider  two  regions  which  differ  only  by 
their  temperature  distribution,  we  can  write  for  the  difference 
in  anelasticity  between  the  two  regions 

e>p[-PvVRrUJ] 

(8) 


where  the  subscripts  E and  W refer  to  the  eastern  and  western 
North  American  models.  Equation  (8)  can  be  solved  for  Tw  to  ob- 
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tain 
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This  expression  can  be  used  to  determine  values  for  E*  and  V* 
which  satisfy  the  temperatures  as  a function  of  depth  for  the  two 
regions,  assuming  that  homogeneous  material  Is  the  same  for  each 
region.  Values  of  AfQ^1),  QoW  and  QqE  are  taken  from  Figures  8 
and  9.  Thermal  models  for  the  eastern  United  States  and  Basin 
and  Range  province  are  given  by  Blackwell  (1971).  These  are  taken 
to  be  representative  of  the  thermal  structures  of  the  uwo  large 
provinces  discussed  above,  for  purposes  of  computation,  although 
the  temperatures  In  the  Basin  and  Range  Province  are  probably 
slightly  higher  than  average  for  western  North  America. 

The  temperature  ratios,  T^^/Tjr  > determined  from  the  values 
of  Blackwell  (1971)  are  shown  as  a dashed  line  In  Figure  10.  The 
value  for  the  activation  volume  has  little  effect  at  crustal 
depths;  consequently  only  an  appropriate  value  for  the  activation 
energy,  E*,  was  required.  Temperature  ratios  for  various  values 
of  E*  appear  In  Figure  10.  It  appears  that  a value  of  1.7 
kcal/mole  could  adequately  explain  the  differences  for  eastern 
and  western  North  America  In  terms  of  temperature  alone.  However, 
this  result  pertains  to  the  relative  values  A(Q_1).  If  the 
same  value  of  E*  Is  applied  to  equation  (7)  for  either  an  eastern 
or  western  North  American  model  alone,  the  values  decrease 
very  rapidly  with  depth  to  values  less  than  80  at  a depth  of  20  km. 
These  are  much  too  small  to  satisfy  the  attenuation  data;  conse- 
quently differences  In  In  the  upper  crust  of  eastern  and 
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western  North  America  cannot  be  explained  solely  by  lateral  vari- 
ations In  temperature  within  a homogeneous  material.  Differences 
In  the  lower  crust  could,  however,  still  be  explained  by  tempera- 
ture dlf  'erences . 
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Figure  Captions 

Figure  1.  Amplitudes  of  Rayleigh  waves  generated  by  Rulison 

(dots)  and  theoretical  radiation  patterns  for  a composite 
explosion  and  double-couple  source . The  four  short  lines 
in  each  diagram  indicate  an  amplitude  of  1.0  x 10“^  cm-sec . 

Figure  2.  Amplitudes  of  Rayleigh  waves  generated  by  Rio  Blanco 
(dots)  and  theoretical  radiation  patterns  for  a composite 
explosion  and  double-couple  source . The  four  short  lines 
in  each  diagram  indicate  an  amplitude  of  2.0  x 10“3  cm-sec. 

Figure  3.  Compresslonal  velocities  (o»0,  shear  velocities  ((3), 
and  densities  (^>)  for  models  of  the  central  United  States 
and  Colorado  Plateau  and  the  Rayleigh  wave  phase  velocity 
curves  which  they  produce. 

Figure  4.  Amplitude  factors,  e"K*  , for  the  models  of  Figure  3. 

Q-L  indicates  the  Qp,  value  in  the  upper  17  km  of  the  crust. 
Q^,  for  the  lower  crust  was  taken  to  be  2000  in  all  cases. 

Figure  5.  Average  values  of  Rayleigh  wave  amplitude  attenuation 
coefficient  ( ^ ) , ground  motion  produced  by  the  explosion 
(We),  relative  strength  of  the  double  couple  component  rela- 
tive to  the  explosion  (F),  and  orientation  of  the  double 
couple  for  the  Rulison  event.  The  vertical  bars  represent 
the  value  of  the  standard  deviation  for  each  determination. 

Figure  6.  The  quantities  described  in  Figure  5 for  the  Rio  Blanco 
event. 

Figure  7.  The  ratio  of  average  amplitudes  recorded  west  and  east 
of  the  sources.  The  solid  and  dashed  lines  result  from  the 
two  models  of  Figure  8. 
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Figure  8.  Values  of  (Q^)  resulting  from  the  solution  of 

equation  (6).  The  solid  line  produces  a very  good  fit  to 
the  data  of  Figure  7 and  the  dashed  line  produces  ampli- 
tude ratio  values  which  differ  from  the  observations  by 
about  30  per  cent  at  short  periods.  Resolving  kernels  for 
several  depths  appear  on  the  right. 

Figure  9.  A modification  of  Q p model  ENA  7 (Mitchell,  1973b) 
for  eastern  North  America  (solid)  and  a model  for  western 
North  America  (dashed)  based  upon  ENA  7 and  the  values  of 
Figure  8. 

Figure  10.  Ratios  of  temperature  needed  to  produce  differences 
between  eastern  and  western  North  American  Q ^ models  for 
various  values  of  activation  energy  (solid  lines),  and  ratios 
derived  from  thermal  models  of  the  eastern  United  States  and 
Basin  and  Range  province  (dashed  line) 
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B.  Work  in  Progress  j 

Work  in  progress  includes  two  studies,  one  of  the  surface 
wave  attenuation  in  Eurasia,  the  ether  of  the  radiation  of  seis- 
mic waves  from  shallow  sources  in  the  Sakhalin  and  Kurile  Is- 
lands  . 
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B.  1.  Surface  Wave  Attenuation  in  Eurasia 
Introduction 

Surface  wave  attenuation  coefficients  determine  the  rate  at 
which  surface  wave  energy  dissipates  as  it  travels  away  from  an 
earthquake  or  an  explosion.  Values  for  the  observed  attenuation 
coefficients  are  determined  in  ideal  cases  by  the  intrinsic  Q 
values  of  the  materials  through  which  the  waves  travel.  In  real 
cases,  however,  attenuation  is  also  produced  by  scattering,  lat- 
eral refraction,  and  mulipathlng  produced  by  non-ideal  source 
conditions,  geologic  complexities,  and  lateral  variation  of  veloc- 
ity and  attenuative  properties  of  the  medium. 

Por  the  reasons  cited  above,  attenuation  data  are  more  dlf- 
fl 'uit  to  obtain  than  velocity  data,  and  consequently  the  amount 
of  attenuation  data  presently  available,  especially  in  certain 
regions,  are  quite  limited.  The  Eurasian  continent  is  one  of 
the  regions  for  which  surface  wave  attenuation  data  are  limited. 
The  only  presently  available  data  for  sources  in  Eurasia  are 
those  of  Tryggvason  (1965)  and  Burton  (197*0-  Both  of  these 
studies,  however,  Include  paths  which  traverse  both  continental 
and  o'eanic  regions;  their  results  consequently  may  not  pertain 
to  true  attenuation  within  Eurasia. 

The  present  work  will  restrict  both  sources  and  stations  to 

the  Eurasian  continent.  This  approach  avoids  complications  which 
arise  from  waves  which  travel  through  ocean  basins  ox*  which  cross 
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continental  boundaries.  The  disadvantage  of  this  approach  is 
the  relatively  small  number  of  stations  available  for  the  deter- 
mination of  attenuation  coefficients. 

We  seek  to  determine  attenuation  coefficient  values  of 
periods  between  10  and  about  100  seconds . The  shorter  periods 
are  especially  affected  by  the  complexitlec  ''cussed  above. 

The  long  path  lengths  from  sources  to  seismograph  stations  in 
Eurasia  combined  with  complexities  along  most  paths  lead  to  rela- 
tively large  standard  deviation  values  for  average  attenuation 
coefficient  values  in  that  region.  Our  approach  will  be  to  pre- 
sent these  average  values  and  their  standard  deviations,  and  then 
to  look  for  regional  departures  from  these  average  values. 


Events 


Data  from  nine  events  are  presently  in  various  stages  of 
analysis.  Table  1 lists  these  events  in  chronological  order. 

The  table  includes  six  earthquakes,  one  underground  nuclear  ex- 
plosion, and  two  atmospheric  nuclear  explosions.  Surface  waves 
from  the  earthquakes  and  underground  explosion  yield  data  in  the 
period  range  from  about  10  to  50  seconds,  and  surface  waves  from 
the  atmospheric  nuclear  events  extend  these  data  to  periods  of 
about  100  seconds  . We  attempted  to  obtain  a wide  distribution 
of  events  throughout  Eurasia.  We  were  limited,  however,  by  the 
absence  of  useful  events  over  much  of  the  region  under  study. 


Stations 

Stations  of  the  World  Wide  Standard  Seismograph  Network 
(WWS3N)  have  provided  most  of  the  data  used  in  the  present  work. 
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Data  from  stations  of  the  USSR  have  been  ordered  through  World 
Data  Center  A but  have  not  arrived  in  the  several  months  since 
the  order. 

The  absence  of  WWSSN  in  the  USSR  and  China  limits  most  of 
the  available  data  from  eastern  Asia  to  paths  through  the  south- 
ern portions  of  the  continent.  The  Novaya  Zemla  nuclear  events 
do,  however,  provide  a few  paths  further  to  the  north  in  eastern 
Asia . 

The  Geophysical  Institute  in  Peking,  China,  has  provided  film 
copies  of  seismograms  produced  by  the  events  of  November  8,  1971 
and  November  18,  1971.  The  records  are  from  Klrnov  type  instru- 
ments with  a flat  response  between  0.3  and  12  seconds.  Data  from 
these  will  therefore  be  limited  to  relatively  short  periods. 

Method 

The  method  for  obtaining  attenuation  coefficient  values  is 
that  of  Mitchell  (197^).  It  assumes  chat  any  event  can  be  con- 
sidered as  an  explosion  plus  a horizontal  double  couple  component 
and  can  be  described  by 

A(uj)  = We  £l  + F sin  2(0  - 0o)  } exp  £ (l) 

where  A(uj  ) is  the  spectral  amplitude,  We  is  the  component  of  mo- 
tion due  to  an  explosion,  P is  the  ratio  of  energy  due  to  the 
double  couple  as  compared  to  an  explosion,  0o  is  the  orientation 
of  the  double  couple,  y'is  the  attenuation  coefficient,  x is  the 
distance  between  the  source  and  receiver,  and  0 is  the  azimuthal 
direction  from  the  source.  The  quantities  We,  F,  and  are 
functions  of  frequency  (to  ) and  have  physical  meaning  only  when 
the  source  is  a composite  of  an  earthquake  and  an  explosion. 
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It  is  found  that  patterns  of  spectral  amplitudes  described 
by  equation  (l)  can  assume  circular,  elliptical,  two-lobed,  or 
four  lobed  shapes,  depending  on  the  value  of  F.  Consequently, 
most  any  radiation  pattern  produced  by  explosions  or  earthquakes 
can  be  described  by  the  equation.  Values  for  We,  F,  ©0,  and  y' 
can  be  determined  from  a least-squares  fit  between  observed  and 
theoretical  radiation  patterns  at  each  period.  Altnough  We,  F, 
and  ©0  will  not  have  physical  meaning  except  in  special  cases, 

Y will  be  the  true  attenuation  coefficient  in  all  cases. 

Results  to  date 

The  Chinese  nuclear  event  of  September  9,  1969,  the  Chinese 
earthquake  of  November  1,  1971.  the  U.S.S.R.  earthquake  of  Octo- 
ber 28,  1971,  the  U.S.S.R  earthquake  of  November  18,  1971  have 
yielded  average  Rayleigh  wave  attenuation  coefficient  values.  In 
each  case  the  values  are  characterized  by  greater  standard  devi- 
ation values  than  those  obtained  from  work  in  North  America.  A 
relatively  small  number  of  available  stations  and  long  paths 
over  geologically  complex  regions  probably  combine  to  produce 
these  standard  deviations. 

The  Chinese  earthquake  produces  low  values  for  the  Rayleigh 
wave  attenuation  coefficient  over  the  entire  period  range  between 
10  and  50  seconds.  The  values  range  from  slightly  negative  to 
about  0.03  x 10"5  km-1  and  have  standard  deviations  of  about 
0.2  x 10~3.  The  Chinese  nuclear  event  produces  values  which  are 
similar  to  those  above  between  10  and  50  seconds,  but  extends 
the  range  of  data  to  about  100  seconds.  The  Rayleigh  wave  atten- 
uation coefficient  values  remain  low  to  periods  of  70  seconds, 
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then  appear  to  increase  to  about  0.3  x 10~3  at  a period  of  100 
seconds.  The  standard  deviation  through  m03  , of  the  period  range 
is  about  0.1  x 10"3. 

The  Russian  earthquakes  exhibit  higher  values,  but  within 
the  error  limits  all  of  the  values  are  similar.  The  values  for 
these  cases  decrease  from  about  0.6  x 10“3  km-1  at  a period  of 
10  seconds  to  aoout  0.2  x 10"3  km"1  at  a period  of  20  seconds 
and  remain  nearly  constant  out  to  a period  of  50  seconds.  The 
standard  deviations  are  about  0.2  x 10~3  for  most  of  the  period 
range . 

Some  of  the  values  listed  above  may  change  if  we  can  supple- 
ment our  present  data  set  with  data  from  stations  in  the  U.S.S.R. 

Regional  variations  in  attenuation  seem  to  occur,  especially 
at  the  shorter  periods.  Consistently  high  amplitudes  are  recorded 
by  stations  in  Fennoscandla  and  are  probably  due  to  relatively 
low  attenuation  across  the  Baltic  snield.  Lower  than  average 
amplitudes  occur  at  Aslan  stations  for  paths  across  the  Himalayas 
and  Tibeten  Plateau.  It  is  uncertain  as  yet  whether  these  low 
values  are  due  to  intrinsically  lower  Q values  in  the  region  or 
to  scattering  effects. 
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B.  2.  Radiation  from  Shallow  Focus  Earthquakes 
in  the  Sakhalin  and  Kurile  Islands. 

A sequence  of  six  earthquakes,  a main  shock  and  five  of 
its  larger  aftershocks,  occurring  in  the  Sakhalin  Islands  in 
September  1971  have  been  selected  for  study.  The  depths  of  fo- 
cus for  these  earthquakes  vary  from  6 km  to  24  km.  Focal  mech- 
anism solutions  have  been  determined  for  each,  and  have  been  found 
to  vary  from  almost  pure  thrust  motion  in  the  main  shock  to  a pro- 
gressively more  dominant  strike  slip  motion  in  the  aftershocks. 

The  dates  of  occurrence,  locations,  and  focal  parameters 
of  these  earthquakes  are  given  in  Table  1 . 

With  the  determination  of  the  focal  mechanisms  as  a base, 
work  is  now  in  progress  in  comparing  the  Ms  vs  mb  values,  particu- 
larly at  Eurasian  stations.  This  will  be  followed  by  a study  of 
the  computed  and  observed  spectral  characteristics  to  examine 
to  what  extent  the  shallow  depth  of  focus  and  thrust  mechanism 
may  affect  determinations  of  Ms. 

Other  examples  of  shallow  thrust  fault  earthquakes  in  the 
Kurile  Islands  (Stauder  and  Mualchin,  1974)  are  available  and 
will  be  used  in  comparison  with,  and  in  supplementing  informa- 
tion from,  the  Sakhalin  Islands  earthquakes. 
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Table  1.  Focal  Parameters  of  Sakhalin  Island  Earthquakes. 
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